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Abstract 

We study the locality of the acceleration temperature in the Unruh effect. To this end, we 
develop a new formalism for the modeling of macroscopic irreversible detectors. In particular, 
the formalism allows for the derivation of the higher-order coherence functions, analogous to 
the ones employed in quantum optics, that encode temporal fluctuations and correlations 
in particle detection. We derive a causal and approximately local-in-time expression for an 
Unruh-Dewitt detector moving in a general path in Minkowski spacetime. Moreover, we derive 
the second-order coherence function for uniformly accelerated Unruh-Dewitt detectors. We 
find that the fluctuations in detection time for a single Unruh-Dewitt detector are thermal. 
However, the correlations in detection-time between two Unruh-Dewitt detectors with the 
same acceleration but separated by a finite distance are not thermal. This result suggests 
that the Unruh effect is fundamentally local, in the sense that the notion of acceleration 
temperature applies only to the properties of local field observables. 



1 Introduction 



A fundamental property of quantum field theories on Minkowski spacetime is that the number 
of particles in a given quantum field state is the same for all inertial observers. However, this 
equivalence does not hold for non-inertial observers, because such observers define particles with 
respect to different normal modes of the field [H [2] . This fact is aptly demonstrated by the Unruh 
effect: for an observer moving with uniform proper acceleration a, the usual quantum field vacuum 
appears as a heat bath at temperature T = The relation between acceleration and temperature, 
characterizing the Unruh effect, has strong analogies to the particle emission from black holes |1] 
and cosmological horizons [5] and for this reason, it constitutes a fundamental ingredient of theories 
advocating a thermodynamic origin of gravity j6]. 
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In quantum field theory, the Unruh effect follows from the fact that the restriction of the 
Minkowski vacuum in one Rindler wedge is shown to satisfy the Kubo-Martin-Schwinger (KMS) 
condition for thermal states [7J. This derivation employs the idealized notion of eternally acceler- 
ated observers and it depends on global properties, such as the existence of a Rindler horizon [8]. 
For this reason, it is important to confirm the physical relevance of acceleration temperature in 
models that involve only local physics [21 El [IHl [HI [121 [13] • The most commonly employed model 
is the Unruh-Dewitt detector, introduced in Ref. [12], which consists of a point-like monopole 
interacting with a quantum field and moving with uniform acceleration in Minkowski spacetime. 

In this paper, we approach the Unruh effect by studying the response of accelerated macro- 
scopic detectors, analogous to ones that have been long employed in quantum optics [T6]. The 
use of macroscopic detector models allows for the consideration of more general observables that 
correspond to the temporal fluctuations and correlations for the recorded particles. In particular, 
we construct the second-order coherence function for accelerated detectors that encodes the fluctu- 
ations of their recorded intensity [17]. We find that the intensity fluctuations recorded by a single 
accelerated detector are thermal, in the sense that the corresponding second-order coherence func- 
tion coincides with that of a static detector in a thermal bath. However, the intensity correlations 
of a pair of Unruh-DeWitt detectors, moving with the same acceleration along different paths are 
not thermal: they differ significantly from the intensity correlations recorded by a pair of static 
detectors in a thermal bath. 

The results above emphasize the local character of acceleration temperature, in the following 
sense. The intensity fluctuations of a single detector correspond to local field observables, while 
the correlation between spatially separated detectors correspond to non-local field observables. 
Therefore, we conclude that the equivalence between acceleration and temperature persists for 
local observables, even when higher-order correlations are taken into account, but it is lost when 
one considers non-local physical observables. 

The consideration of macroscopic detectors is crucial for the derivation of the correlations and 
fluctuations of the recorded particles, and it is the key point of departure of this paper from 
past works on the topic, where detectors are treated as microscopic systems. Microscopic systems 
are meaningful probes of a quantum field, and as such they suffice to ascertain the existence of 
acceleration temperature, but they do not conform to the conventional concept of a detector [H]. 
An actual detection event involves an irreversible amplification of microscopic processes, a definite 
macroscopic record and it is well localized in time. 

Using only microscopic probes, it is impossible to associate a definite instant of time (or time 
interval) with a detection event. An elementary example that demonstrates this point is the fol- 
lowing. Let us consider a single (isolated) atom interacting with a quantum field through a dipole 
coupling. In absence of a measurement, i.e., an interaction with a macroscopic measuring device 
that leaves a definite macroscopic record, we cannot tell when during its history, the atom absorbed 
a photon. A definite record requires that (at least part of the) detecting system be macroscopic. 
One may consider for example a detector consisting of an array of atoms, together with a macro- 
scopic photo-detector. When an atom absorbs a photon, it becomes excited. However, the excited 
state has a finite lifetime F"^; the atom then re-emits a photon which leaves a macroscopic record 
at the photo-detector. Thus, the moment of detection can be localized in time with an uncertainty 
of order F~^. In absence of the macroscopic photo-detector, the instant of detection cannot be 
specified. 

In existing treatments, the response of a microscopic probe coupled to a scalar field is de- 
termined from the solution of Schrodinger's equation for the total system. To first order in 
perturbation theory, the transition amplitude for a uniformly accelerated Unruh-Dewitt detec- 
tor corresponds to a Planckian spectrum [15] . Higher-order effects, such as dissipation and noise, 
destroy the thermal behavior for proper times much larger than the decay time scale of the detector 
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|18] . This problem does not arise if one assumes that the detector- field interaction is switched on 
only for a finite time interval [THl EHl EI] of duration T much smaller than the decay time scale. 
In this case, the transition amplitude refers to the state of the detector after the interaction has 
been switched-off: it does not define the detection probability during the period of interaction. To 
first order in perturbation theory, the transition amplitude is proportional to the time T, hence, it 
is possible to infer an expression for the detector's excitation rate at intermediate times. With a 
suitable regularization, one obtains a causal expression for the excitation rate of a detector moving 
along general trajectories in Minkowski spacetime [22t |23 | l2l]. 

However, a purely unitary evolution of a microscopic system cannot account for the irreversible 
macroscopic amplification corresponding to a detection event and the creation of an associated 
record. This implies, in particular, that correlations between detection events cannot be deter- 
mined. For example, the joint probability that a particle is detected at time ti and another particle 
is detected at time ^2 cannot be constructed. In order to construct such joint probabilities, it is 
necessary to incorporate the information obtained by the first detection into the quantum state — 
as for example, when one considers sequential measurements p5 | [33l [27] . 

Multi-time correlations are observable in quantum optics [17], where the statistical analysis 
of photon-detection times provides significant information about the quantum field state. The 
photon statistics may be Poissonian, sub-Poissonian, or super-Poissonian, and the photons may 
either bunch or anti-bunch. (See, Ref. [17] for a definition of these terms, and Sec. 3.2.2 for their 
adaptation in the present context.) The measurement of such correlations requires macroscopic 
detectors with well-defined records for the instants of detection. These properties must also be 
refiected in the theoretical description of the detection process. 

To this end, we employ techniques from the quantum theory of measurement [23 ESI EH] 
in order to construct a general model for macroscopic particle detectors, moving along general 
trajectories in Minkowski spacetime. In particular, we adapt a method that was developed in 
Refs. [301 EI] foi' the study of time-of-arrival and related problems in quantum theory. This 
method contains some ideas and techniques from the consistent /decoherent histories approach to 
quantum mechanics [321 [2Sl ESI EH ES], and it is related to the Davies-Srinivas photodetection 
theory [36]. The most significant advantage of this method is that it leads to the definition of 
detection probabilities where the time of detection t is treated as a genuine quantum variable, and 
not as an external parameter. In particular, the detection probabilities are densities with respect to 
t and they define a Positive-Operator- Valued-Measure (POVM). The treatment of detection time 
as a quantum observable is essential for an unambiguous and probabilistically sound definition of 
the fiuctuations and correlations with respect to the time of detection. 

Using this method, we obtain general expressions for the detector's response, which are causal 
and approximately local in time. We also derive expressions for the multi-time joint probabilities 
of particle detection. The correlations for a pair of uniformly accelerated detectors are thermal if 
the detectors are co- incident, but non-thermal if they are separated by a finite distance d. 

The structure of the paper is the following. In Sec. 2, we develop a general formalism that 
determines the intensity and the intensity fiuctuations of irreversible particle detectors. In Sec. 
3, we apply this formalism to detectors with a dipole coupling to a massless scalar field (Unruh- 
DeWitt) detectors. We derive expressions for the detection probabilities and correlations for such 
detectors moving along general trajectories in Minkowski spacetime, and we then specialize to the 
case of constant acceleration. In Sec. 4, we summarize our results, emphasizing the local character 
of acceleration temperature. 
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2 Probabilities for an irreversible particle detector 



In this section, we develop a general formalism that determines the joint probabilities for a sequence 
of quantum events, with respect to the times occurring at different time instants ti,t2, ■ ■ ■ ,tn- 
The novelty in this construction is that the times ti,t2, ■ ■ ■ ,tn are treated as genuine quantum 
observables. The probabilities are densities with respect to (ti,t2, ■ ■ ■ ,tn) and they correspond to 
a Positive-Operator- Valued-Measure. These probabilities are then employed in order to derive the 
intensity and the intensity fluctuations of irreversible, macroscopic particle detectors. 

The main requirement in this construction is the existence a coarse-graining time-scale a, such 
that all temporal interferences at time-scales larger than a are suppressed. In the language of 
decoherent histories, this means that histories corresponding to a transition taking place at some 
time interval [ti,tf] are decoherent, if \tf — ti\ » a. This property is expected to hold for 
macroscopic systems, or systems that possess a macroscopic component [32l |26| |3^. 

2.1 Probability amplitudes for the time of detection 

The first step in our construction is the derivation of a general expression for the single-time 
probability for a quantum event taking place at time t. The event time t is treated as a quantum 
observable, and the probabilities are densities with respect to t. The more general case of multi- 
time probabilities associated to a sequence of quantum events is treated in Sec. 2.3. 

First, we derive an expression for the probability amplitude that, given an initial state iV'o)) a 
transition occurs at some instant in the time interval and a recorded value A is obtained for 

some observable. 

Let 7/ be the Hilbert space of a quantum system. Since we are interested in describing mea- 
surements, % includes the degrees of freedom of a macroscopic apparatus in addition to the ones 
of a microscopic system. In order to describe a quantum event, like a particle detection, we assume 
that % splits into two subspaces: 1-i = 71+ The subspace describes the accessible states 

of the system given that a specific event is realized; the subspace "H- is the complement of "H^. For 
example, if the quantum event under consideration is a detection of a particle by a macroscopic 
measuring device, describes the possible states of the system provided that the detection has 
occured and "H. describes the possible states of the system if the detection has not occurred. We 
denote the projection operator onto "H^ as P and the projector onto Ti^ as Q := 1 — P. 

Once the transition has taken place, it is possible to measure the values of various observables 
through their correlation to a pointer variable. We denote by Pa projection operators corresponding 
to different values A of an observable that can be measured only if the quantum event under 
consideration has occurred. For example, when considering a particle detection, Px may correspond 
to possible values of a pointer variable in the apparatus that is correlated to some microscopic 
property of the system: position, momentum, energy, and so on. The set of projectors Pa is 
exclusive {P\Px> = 0, if A 7^ A'). It is also exhaustive given that the event under consideration has 
occurred; i.e., ^a-^a = P- 

We assume that the system is initially (at time t = 0) prepared in a state \ipo) G and that 
the dynamics is governed by the self-adjoint Hamiltonian operator H. 

We proceed to derive an expression for the amplitude I'?/'; A, [^1,^2]), that a transition occurs 
at a time interval [^1,^2] and a recorded value A is obtained. We first consider the case that the 
relevant time interval is small, i.e., we set ti = t and t2 = t + 6t, and we keep only terms of leading 
to 6t. Since the transition takes place within the interval [t,t + 6t], at times prior to t the state 
lay within T-L-. This is taken into account by evolving the initial state \tpo) with the restriction of 
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the propagator into that is, with the operator 

St = hm (Qe-^^*/^Q)^. (1) 

By assumption, the transition took place some time within the time interval [t, t + 6t] and at 
the end the outcome A is recorded. This means that in the time-interval [t, t + 6t\ the amplitude 
transforms under the full unitary operator for time evolution e"*^"^* ~ 1 — i6tH. At time t + 6t 
the event corresponding to Px is recorded, so the amplitude is transformed by the action of P\ 
(or of VPx, if -Pa is not a projector). For times greater than t + 6t, there is no constraint, so the 
amplitude evolves unitarily until some final moment T. 

At the limit of small 6t, the successive operations above yield 

lijo; A, [t,t + St]) = -i5t e-'^^^-'^PxHSt\iJo). (2) 

We emphasize here the important physical distinction on the role of time that is manifested 
in the derivation of the amplitude Eq. ([2]). The time of detection t does not coincide with the 
evolution parameter of Schrodinger' s equation. Instead, it is a physical observable that corresponds 
to the moment that an event has taken place. The construction of the amplitude {ip; A, [t,t + 6t]) 
above highlights this distinction: the detection time t is distinct from the time T at which the 
amplitude is evaluated. 

The amplitude {ipo] A, [t, t + 6t]) is proportional to 6t, hence it defines a density with respect to 
time: |V^o;A,t) := limst^o j^lipo] X,[t,t + 6t]). From Eq. ([2]) 

|V^o;A,t) = -te-'^^^-'^PxHSt\i^o) = -te~'^^Ci\,tMo), (3) 
where the class operator C{X,t) is defined as 

CiX,t) = e'^'PxHSt. (4) 

Since the amplitude \ipo', A, t) is a density with respect to the time of transition t, its integration 
with respect to t is well-defined. Hence, the total amplitude that the transition occurred at some 
time in the interval [^1,^2] is 

-iHT 



|7/>;A,[ti,t2]) = -^e-^^' / dtC{X,t)\^,). (5) 



We note that if [P, H] = 0, i.e., if the Hamiltonian evolution preserves the subspaces H^, then 
\ipQ; A, t) = 0. For a Hamiltonian of the form H = Hq + Hj, where [Hq, P] = 0, and Hi a perturbing 
interaction, we obtain to leading order in the perturbation 

^(A, t) = e'^''PxHie-'^°K (6) 
Eq. ([6]) also holds for a time-dependent interaction Hamiltonian Hj. 

2.2 Probabilities for detection time 

The amplitude Eq. squared defines the probability p(A, [ti, ^2]) that at some time in the interval 
[^1,^2] a detection with outcome A occurred 

p(A, [ti,t2]) := A, [ti,t2]\^; A, [t^^t^]) = T dt T dt' Tr{e'^^'-'"^ PxHSlp^S,HPx), (7) 
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where Po = |V'o)(V'o|• 
However, the expression p(A, [^1,^2]) does not correspond to a well-defined probability mea- 
sure, because it fails to satisfy the additivity condition. To see this, consider the probability 
corresponding to an interval [^1,^3] = [^1,^2] U [^2,^3]- This equals 



p{\, [h, h]) = piX, [h, t2]) + piX, [t2, h]) + 2Re 



t2 



dt ' dt'Tr (^C{X, t)poC\X, t') 



Hence, the additivity condition p(A, [^1,^3]) = p(A, [^1,^2]) +p(A, [^2,^3]), necessary for a consistent 
definition of a probability measure, fails, unless 



2Re 



^ dt r dt'Tr (C{X, t)poC\X, t'] 

tl Jt2 ^ 



(9) 



In the consistent / decoherent histories framework, Eq. dH]) is referred to as the consistency condition. 
It is the minimal condition necessary for the association of a consistent probability measure in 
histories. It appears in the present framework, because we consider properties of the system at 
different moments of time, that is, histories. 

Eq. ([9]) does not hold for generic choices of ti, ^2 and ^3. However, in a macroscopic system (or 
in a system with a macroscopic component) one expects that Eq. ^ holds with a good degree of 
approximation, given a sufficient degree of coarse-graining. This means that there exists a time- 
scale 0", such that the non-additive terms in Eq. ([8]) are strongly suppressed if \t2 — ti\ » a and 
1^3 — ^2! >> o"- Then, Eq. ([7]) does define a probability measure when restricted to intervals of size 
larger than a. 

In principle, the coarse-graining time-scale a can be estimated from the knowledge of the 
Hamiltonian and the initial state of the system. Its value depends strongly on the detailed physics 
of the measurement scheme, and constitutes an inherent feature of a measurement process. Using 
a semi-classical language, a can be said to represent the time necessary for a microscopic process 
to amplify and leave a definite macroscopic record. Thus, a places an absolute lower bound in 
the temporal resolution of each measurement scheme: it is impossible to distinguish events if their 
separation in time is of order a or less. 

Temporal smearing. Assuming a finite coarse-graining time-scale cr, such that Eq. ([9]) is 
approximately valid for \t2 — ti| >> o" and 1^3 — ^2! >> cr, Eq- © provides a consistent definition 
of a probability measure. However, it is more convenient to proceed by smearing the amplitudes 
Eq. ([3]) at a time-scale of order a rather than employing probabilities for sharply defined time- 
intervals, as in Eq. ([7]). Smearing allows us to express the probabilities in terms of densities with 
respect to a continuous variable. 

To this end, we introduce a family of functions /cr(s), localized around s = with width a, and 
normalized so that \im.T-^Q f^{s) = S{s). For example, one may employ the Gaussians 

Us) = -L=e-^. (10) 



V2 



The Gaussians Eq. (ITUl) satisfy the following equalities. 

Us)W) = U (^) Us - .') (11) 

^Ut-s)Ut-s') = Ut - '^)9As - s'), (12) 
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where 



g,{s) = exp[-sV(8a2)]. (13) 

Eqs. ( ITTIfT2|) are satisfied approximately for otlier choices for the smearing functions with an error 
of order cr^. 

Using the functions /o-, we define the smeared amphtude iV'o; ^i i)a that is locahzed around the 
time t with width cr, as 

|V'o;A,t).:= j ciV/a(s-t)|^o;A,s) = a(A,t)|^o), (14) 

where 

a(A,t) := j ds^U(t-s)C{Ks). (15) 

The square amphtudes 

p(A, t) = ^(^o; A, tlV^o; A, t), = Tr [cl{\, t)p^C„{\, t)] (16) 

provide a well-defined probability measure: they are of the form Tr [/)on(A, t)], where 

Ii{\t) = C,{\t)Cl{\t) (17) 

is a density with respect to both variables A and t. 
The positive operator 

dt j dXIiri^t), (18) 

corresponds to the alternative M that no detection took place in the time interval [0, oo). nT-(A^) 
together with the positive operators Eq. ( 1T7|) define a POVM on ([0, oo) x fi) U {A/"}, where VL is 
the space of possible values of A. The POVM Eq. ( ITTl) determines the probability density that a 
transition took place at time t, and that the outcome A for the value of an observable has been 
recorded. 

2.3 Multi-time joint probabilities 

Next, we generalize the results above, by constructing the joint probability for transitions that 
occur at different moments of time. 

First, we consider the case of two detection events at times ti and t2 (^i < t2). Again, we 
denote the system's Hilbert space as "H. Let the transition at time ti {i = 1,2) correspond to a 
splitting of "H into two orthogonal subspaces associated to the projectors -Pj+ and Pj_ = 1 — Pj+. 
At each time ti, the values Aj are being recorded; these values correspond to an exclusive set of 
projectors P\., such that J2xi P>^i = 

The probability amplitudes l^/'o; Ai, [ti, ti + 5ti]; A2, [^2, ^2 + 5^2]) that correspond to the recording 
of values Aj at times ti are obtained through a similar procedure to the one leading to Eq. ([2]): 

iVo; Ai, [ti, ti + 5ti]- A2, [t2, t2 + 5t2]) = -6ti5t2e-'''^CiXi, ti; A2, t2)|^o), (19) 

where 

C'(Ai,ti; A2,t2) := e'^'-P,X-tAiHSll'\ (20) 
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('12') 

In Eq. ([20D, is the restricted propagator corresponding to no transitions (either with 

respect to Pi+ or with respect to -P2+), that is, 



Sl''^ = hm 



p -iHt/n p 
-'12-6 -ri2- 



(21) 



where pu- is the projector corresponding to the logical disjunction of Pi- and P2-. For [A-? A-] = 
0, which is the case of interest in this paper (it corresponds to independent apparatuses), -P12- = 
P1-P2-. A similar expression defines S"^ , the restricted propagator corresponding to a transition 
with respect to -Pi+, but not with respect to P2+- 

The explicit form of the restricted propagators is not important when the Hamiltonian H is of 
the form H = Ho + Hj, where [Hq, Pj+] = and Hj is a small perturbation. In this case, the class 
operator C*(Ai,ti; A2,t2) of Eq. (I2U]) becomes 

C'(Ai,ti; A2,t2) = e'^"''Px,Hie-'^''^''-''^Px,e-'^'''\ (22) 

to leading order in perturbation theory. 

To obtain the associated probability density, we smear the amplitude f|T9|) at each moment of 
time ti and t2, using smearing functions fa-^ and fa-2, respectively. In general, the coarse-graining 
time-scales (jj may be different for each transition. Following the same procedure as in Sec. 2.2, 
we obtain 

p(Al,t2; A2,t2) = y dSids[dS2dS2^y faiih - Si)/^i(ti - S[)fa2it2 " S2)fa2ih " ^2) 

X {tPolC^Xu si; A2, S2)C{Xu s[; A2, 4)lV'o). (23) 

The procedure above can be applied systematically in order to obtain the n-time joint prob- 
abilities p(Ai, ti, A2, t2; . . Xn, in) that correspond to measurements of observables Ai, A2, . . . , A„ 
at times ^1,^2, ■ ■ ■ ,tn respectively (ti < ^2 < • • • < tn)- Here, we give the final expression for a 
Hamiltonian Hq + Hj, where [Hq, Pj+] = and Hj is a small perturbation, to leading order in 
perturbation theory 

/n 
Yl dSids'i \/ faAti - Si)fai{ti - S[) 
i=l 

n n 



j=i k=i 



2.4 A general model for particle detectors 

So far, we have constructed expressions for the single-time and multi-time probabilities for a 
generic quantum system. In what follows, we study the case of an apparatus that responds to 
the excitations of a microscopic system and records their energy. We obtain expressions for the 
intensity of the absorbed energy as a function of time and for the intensity fluctuations. 



2.4.1 Intensity at a single moment of time 

We consider a Hilbert space T-Ldet, describing the degrees of freedom of a macroscopic detector, 
and a Hilbert space Hq, describing a microscopic quantum system interacting with the detector. 
The Hilbert space of the combined system is then = Hdet ® Tig. We denote the Hamiltonian 
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of the detector degrees of freedom as Hd and the Hamiltonian of the measured system as Hq. A 
general form of the interaction Hamiltonian is A" ® -Bf , where A"- are operators on T-Ldet and 
131 (possibly time-dependent) operators on "Hg. 

Our aim is to model a detector that records the energy corresponding to the excitations of the 
quantum system. The relevant transitions correspond to changes in the detector's energy, as it 
absorbs an excitation (particle) of the microscopic system. We assume that the initially prepared 
state of the detector has minimum energy Eq, conveniently taken equal to zero, and that the 
excited states have energies E > Eq. We denote the projectors onto the constant energy subspaces 
as Pe- Thus, the projectors P± on "H corresponding to the transitions are P+ = (^e>Eo -^^) ® ^ 

and P- = Po ® i, where Pe = Peq- 

The class operators Eq. for this system are 

CiE, t) = J2 (^i^^"(i)) ® B'^it), (25) 

a 

where A'^{t) = ^^^dt ^a^-iH^t ^-^^ B'^[t) = Q^^qtBle'^^i^ are the Heisenberg-picture evolution of A" 
and B1 respectively. 

Physically reasonable detectors are characterized by an energy gap, i.e, there exists a minimal 
value Emin > Eq in the energies of the excited states. If the detector's temperature is much smaller 
than Emin, then a thermal state for the detector' degrees of freedom is well approximated by the 
density matrix = Pq/TtPq. 

We consider an initial state of the total system p'^^^ ® po, where po is the initial state of the 
microscopic system. Eq. flTBl) then gives 

p{E,t) = J]a»'(i?) j ds j ds'^JUt - s)Ut - s')e-'^^''-'^Trn, [poB'^is')B''{s)] , (26) 
where 

Trn,,, (PoA-PeA") 
^ab ^ V ^ J__ (27) 

The intensity / recorded by the detector at time t (energy absorbed per unit time) is then 
obtained from Eq. ( 126|) 



/(t) = j dEEp{E,t). {2t 



2.4.2 Intensity fluctuations 

Next, we consider the intensity fluctuations measured by one or more detectors. To this end, we 
assume that the detector, being macroscopic, consists of a large number of independent subsystems, 
each of which being able to record a microscopic excitation. The intensity refers to the totality of 
detection events, irrespective of the subsystem where they have been recorded. Hence, in order to 
construct the n-time intensity correlation function, we must consider the response of n independent 
detectors, identical to the ones considered in Sec. 2.3. 

We first work for n = 2; the generalization to higher n is straightforward. We consider two 
copies of a detector's Hilbert space Tidet, so that the Hilbert space of the total system is Tideti ® 
'Hdet2®'Hq. The projector P+i corresponding to the transition recorded by detector 1 is 'Y1,e>Eo Pe® 
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1®! and the projector P+2 corresponding to the transition recorded by detector 2 is ^®^e>Eo 
1. 

We assume that there is no direct interaction between the two detectors. The general form of 
the interaction Hamiltonian is then 

= ^ (^i" ® i ® 5," + i ® ® B^^ . (29) 

a 

Substituting into Eq. fl2^ we obtain the following expression for the class operator C{Ei,ti, -^2,^2) 



a,b 



(30) 



In the derivation of Eq. fl3Ql) . we assumed that the first detection (at time ti) was made by the 
detector 1 and the second detection (at time ^2) by the detector 2. If the time-ordering of the two 
detections is not observable, we must also take into account a contribution to the amplitude with 
the reverse time-ordering for the detections. Then, the relevant class operator is 



a.b 



(31) 



where T stands for time-ordering. 

Thus, from Eq. f l2^ we find the probability density p{Ei,ti; -^2,^2) for a detection with energy 
El at time ti and a detection with energy E2 at time t2 



p{Ei,ti,E2,t2) = Yl a'''\Ei)a^'\E2) ! dsids'^ds^ds'^e''^^^''^-'^ 

a,b,a'b' 



X 



VUlitl - Sl)/<Ti(tl - s[)U^{t2 - S2)fa2{t2 - S'2, 



p,A{B'^\s\)&\s'^)]T{B\s2)B'^{s 



where A stands for the anti-time-ordered product. 

The two-time intensity correlation function is defined as 

(/(ti)/(t2)) := j dEidE2EiEip{EiM\E2M) 
We similarly derive the n-time detection probability density 

n 

p{Ei,h-E2,t2;...;En,tn) = Y Yl n(""'"^(^*) 

ai,...,a„ a[,...,a'„ i=l 



(32) 



(33) 



dsids[e-'^^^<-'^^ ^/f^M^7i)fAt~^ 



X Trn, [poA(^B<{s[)...B<{s'j)r[B^-{sr,) 
and the n-time intensity correlation function 

(/(tl) . . . J(t„)) = J dEi... dEnEi . . . EnPiEi, h; E2, ^2; • • • ; ^n, tn). 



,(34) 



(35) 
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3 The Unruh-DeWitt detector 



In Sec. 2, we derived the equations for the joint probabihties and intensity correlations corre- 
sponding to general macroscopic detectors interacting with a quantum system. In this section, we 
specialize to detectors of the Unruh-DeWitt type that move along general trajectories in Minkowski 
spacetime and interact with a massless scalar field. 



3.1 Detection probability for Unruh-DeWitt detectors 
3.1.1 Derivation of a general expression 

An Unruh-DeWitt detector is an ideal particle detector that interacts with a quantum field via a 
dipole coupling. One usually considers such detectors moving along non-inertial paths in Minkowski 
spacetime. In this context, a crucial assumption is that the physics at the detector's rest frame 
is independent of the path followed by the detector [14]. Hence, the evolution operator for an 
Unruh-DeWitt detector is e"*^**^, where r is the proper time along the detector's path and Hd is 
the Hamiltonian for a stationary detector. Thus the path-dependence enters into the propagator 
only through the proper time. In particular, supposing that the detector follows a timelike path 
x^{t) in Minkowski spacetime, the equation x^{t) = t can be solved in order to determine the 
proper time r(t) as a function of the inertial time coordinate t of Minkowski spacetime. Then the 
evolution operator for the detector with respect to t is e^^^'^'^^^\ 

We assume that the detector is coupled to a massless scalar field 0. The scalar field Hamiltonian 
is Hff) = f d?x (Itt^ + i(V0)^). The evolution operator for the field with respect to a Minkowskian 

inertial time coordinate t is e"*"^"^*, where t is a Minkowski inertial time parameter. The Unruh- 
DeWitt detector interacts with the field through a dipole coupling. This means that the (time- 
dependent) interaction Hamiltonian is of the form Hi = m®0[x(r)], where 0(x) are the Heiseberg- 
picture fields, and m is an operator on the detector's Hilbert space analogous to the magnetic 
moment in electrodynamics. 

The detection probability for an Unruh-DeWitt detector follows from Eq. (1261) . with the sub- 
stitutions A"" — )■ m, B"'{t) — )■ 0[a;(r)]. Using the detector's proper time r as a time parameter, we 
obtain 

piE,T) = a{E) JdsJ ds'^JUt - s)Ut - s')e-''^(^'"'^)A+[a;(s'),a;(s)], (36) 
where a{E) = Tr-y^^{pErhPom) /Tr-^^pQ, and 

A+(a;, x') = Tr U(x)0(x')pol ^'^^^ 
is the positive-frequency Wightman function. When the field is in the vacuum state, 

^^(^'^') = 7-^773^ -fv h-,. „A.i > (38) 



47r2[(a;0-x0' -ie) 



X — X' 



where e > is the usual regularization parameter. 

Eq. fl5Bl) simplifies using Eq. flTT]) and then setting /^-(r — ^i^) ~ 5{t — ^^)- This is a good 
approximation when the system is monitored at timescales much larger than a. Then, 

p{E,T) = a{E) I dyg^{y)e-^^yA^[x{T + ^),x{t - ^)]. (39) 
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Eq. f l5^ determines the detector's response for motion along a general spacetime path. It is similar 
to the standard expression obtained through first-order perturbation theory with an adiabatic 
switching on and off of the coupling [22| |23l |2l]. The difference lies on the factor ga{y)- This term 
guarantees that the detector's response to changes in its state of motion is causal and approximately 
local in time at macroscopic scales of observation. To see this, we note that ga{y) truncates 
contributions to the detection probability from all instants s and s' such that |s — s'| is substantially 
larger than a. Thus, at each moment of time r, the detector's response is determined solely from 
properties of the path at times around r with a width of order a. In particular, properties of 
the path at the asymptotic past (or future) do not affect the detector's response at time r. The 
response is determined solely by the properties of the path at time r, within the accuracy allowed 
by the detector's temporal resolution. 

At the limit a oo, gcr{y) 1, and the probability p{E, t) is non-local in time and non-causal. 
According to the discussion in Sec. 2.2, this limit corresponds to a physical system for which the 
detection event is not be localized in time. A finite value of a is indicative of a detector with 
many degrees of freedom, and thus the records of particle detection are coarse-grained observables. 
Coarse-graining typically implies noise that may deform the 'signal' due to the detector's state of 
motion. We will see that this is indeed the case for uniform accelerated detectors, where a finite 
value of (J implies a deformation of the Planckian spectrum. 



3.1.2 Uniformly accelerated detectors 



Next, we apply Eq. fl36|) . in order to obtain the response of a uniformly accelerated detector. We 
consider a spacetime trajectory of the form x^{t) = sinh(ar)/a, a;"'^(r) = (cos(ar) — l)/a,x^(r) = 



0. Then, 



piE, 



167r2 J smh\ay -il)/2\ 



(40) 



where a time-dependent factor has been absorbed into e. 

The integrand in Eq. fHU]) is peaked around y = 0. The hyperbolic sinus in the denominator 
has a width of order 2/a around y = 0, and the function g^j has a peak of width a around y = 0. 
If aa >> 1 and E/a is not much smaller than unity, we obtain p{E, r) to leading order in l/{aa) 
by setting g^ = 1. We find 



p{E,t) 



dye 



-iEy 



a{E)a^ 
167r2 J sinh^[a(?/ - ze)/2] 



a{E) E 
2tt e~ - 1 



(41) 



a 

27r" 



corresponding to a Planckian spectrum for the distribution of energy, with temperature T 

To compute corrections due to the finite value of l/{aa), we use Eq. ( fT3l) for g^^ and express 
Eq. as 



p{E,t) 

where r] = (cra)~^ and z = ay. 
We evaluate the integral Eq. 



a{E)a 
' 167r2 



dz- 



-■qz^ /8~iEz/a 

' sinh^ [z -ii)/ 2]' 



(42) 



p{E,r) 



a{E) E 



27r 



2tvE 

e a 



as a power series in rj, to obtain 

2nE 



1 



TC 



e— + 1 



, 2wE 

[e a 




(43) 



Hence, provided that a/E is not much larger than unity, the effect of a finite value for l/(cra) 
is a small deviation from the Planckian spectrum. However, when the condition aa >> 1 is not 
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satisfied, the response function takes an entirely different form and lias no similarity to thermal 
radiation. We note that the condition aa » 1 for the identification of the Unruh temperature is 
mirrored in the study of the Unruh effect for observers with finite lifetimes tq [371138] . The physical 
context in the above references is rather different, but they confirm the point that the derivation 
of the Unruh temperature requires that the time-scale 1/a associated to the acceleration is much 
smaller than any intrinsic time-scales of the detecting system. 

We conclude that, in a macroscopic detector, Unruh's relation between acceleration and tem- 
perature requires a separation of timescales. If the resolution timescale is of the order of 1/a, the 
quantum coherence of the microscopic processes related to particle detection is lost. However, 
even if the resolution timescale a is much larger than 1/a, it may still be much smaller than the 
macroscopic timescales characterizing the detector's trajectory. For example, consider a path char- 
acterized by a stage of uniform acceleration of proper-time duration equal to At » a. During 
this period, the detector's response will be Planckian, modulo transient effects of order a/ At and 
spectrum deformations of order (cra)"^. Thus the local-in-time thermal response of an accelerated 
detector arises in this calculation as a consequence of the separation of three time-scales: the mi- 
croscopic scale 1/a, the resolution scale a and the macroscopic scale At characterizing variations 
in the detector's acceleration. 

We next compare the response of the accelerated detector in the Minkowski vacuum, with 
a response of static detector (x" = t,x^ — 0) in a thermal bath at temperature T = The 
corresponding positive-frequency Wightman function for a massless quantum field is 

/^ ^ I sinh(27rr//3) 

' ' ' ^ 87r/3rsinh[7r(t-f-ze- |x-x'|)//3]sinh[7r(t-f -ze+ |x-x'|)//3] 

oo/ 

J2 S{t-t' + inP), (44) 



oo/ 



An 

n=—oo 



where the prime denotes that the n = term is excluded from the summation. For x = x' 
Eq. (jH]) becomes 



Substituting Eq. (HSl) into Eq. (!39|) . we find that at the limit aa » 1 {go- — 1, the delta 
functions do not contribute to the total integral because they vanish along the real axis. Hence, 
the response function p{E, r) coincides, as expected, with that of an accelerating detector for 
T = This equivalence may be valid for state other than thermal. Any Wightman function 
with poles in the negative imaginary plane at t — t' = —iriK, where n = 1, 2, . . . and k a positive 
constant gives rise to a Planckian spectrum (for na » 1). 



3.1.3 Detector's response for other paths 

The response function Eq. fl5^ also applies to other choices of spacetime paths. We examine some 
cases: 

1. For paths such that A'^[x(s), x{s')] is a function only of s — s', the detection rate is constant. 
This result also follows from first-order perturbation theory in a microscopic detector |14j . 

2. For paths with time-dependent linear acceleration a(r), such that the acceleration varies 
at a scale much larger than a {a /a « 1/cr), the detection rate corresponds to a time-dependent 



13 



Planckian spaectrum (as long as a(r)cr >> 1), i.e., 

p{E,t)^^^ . (46) 

Thus the detector comes to equilibrium at the instantaneous values of the Unruh temperature 
T(r) = ^. 

3. For paths of the form {x^{t),x^{t),0,0) that correspond to motion along a single axis, we 
define u(r) = x°(r) + x^{t) and v{r) = x°(r) — x^{r). Clearly iiv = —1. Eq. (j39|l becomes 

(F \ = -'^!d gMe"'''' 

'^^ 47r2 y ^[«(r + |//2)-«(r-i//2)-ze][t;(r + y/2)-t;(r-y/2)-^e]- ^ ^ 

The function (70- (l/) effectively restricts integration to values of y of order of a or smaller. Thus, 
we can Taylor-expand the terms in the denominator of Eq. (147|) as u{t + y/2) — u{t — y/2) = 

yu{T) + XlfcLi 2^(2fc+i")! ^^^]' ^^"^ truncate the series to some order km if M(2'^'"+'^)(r)cr^'' becomes 
significantly smaller than the previous terms in the expansion. Similarly, we can truncate the 
series for u{t + y/2) — u{t — y/2) at some order k'^, so that the denominator in Eq. f H7|) is 
approximated by a polynomial of order k^ + k'^. We then obtain the dominant contribution to 
p{E, t) by letting g„ ^ 1, whence the integral corresponds to the sum of the residues for the poles 
lying in the negative imaginary plane. For example, if we truncate the series at k^ — ^In — ^ ' the 
integrand in Eq. ( H7|) has a double pole at y = and also poles at y"^ = —2Au/'u and y"^ = —2Av/'v. 
If < ^ v/'v > 0, there are two single poles in the negative imaginary plane, whence 



47r iiv — iju 

In contrast, if ii/'u^v/'v < 0, all poles are real- valued and thus p{E,t) vanishes. 

3.2 Intensity fluctuations for the Unruh-DeWitt detector 
3.2.1 The general case 

Next, we construct the joint detection probability p{Ei, ti, E2, T2) for a pair of identical detectors 
that move along different spacetime trajectories, Xi{ti) and X2(t2). Eq. ( !32l) applies to this 
case. We substitute A" — )■ m, i?"(t) — >■ 0[xi(ri)], B^{t) — )• 0[x2(t2)], and we use the proper-time 
parameters ri and T2, corresponding to Minkowskian times ti and ^2 via Xi{ti) = ti and x^(r2) = t2. 
We obtain 



p(Ei,Ti; E2,T2) = a{Ei)a{E2) J dsids[ds2ds'2e^ 



iEi{s[-si)-iE2{s'2S2) 



X V/a(tl - Sl)faiti - S[)fait2 - S2) f a{t2 - 



X Tr 



Hi 



A U[xi(s;)]0[x2(4)] r 0[x2(s2)]0[xi(si)] po , (49) 



where time- and anti-time-ordering refers to the Minkowski time parameters that are associated to 
the proper times si,S2 and s'^,S2 respectively. When po is the Minkowski vacuum state, Eq. ( H9|l 
becomes 



p(Ei,Ti; E2,T2) = a(Ei)a{E2) J dsids[ds2ds2e~ 



■iEi{s[-si)-iE2{s'2-S2) 



X\/ f<j{tl - Sl)fa{tl - S[)f„{t2 - S2)fa{t2 " 4) ( [0:2 (4) > (s'^ )] Ai;' [X2 (^2) , (si )] 

+A+[a;2(4),a:2(s2)]A+[xi(s;),a;i(si)]+A+[a;2(4),^i(^i)]A+[a;i(s;), 0:2(52)]) , (50) 
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where Ap{x, x') = —1/[Att^{x — x')^] — — x')"^] is the Feynman propagator for the scalar field. 

The second product of Green's function in Eq. (11^ gives rise to p{Ei, Ti)p{E2, . Thus, 

p{Ei, n; E2, ra) = p{Ei, n)p{E2, T2) + G{Ei, T2; E2, T2), (51) 

where the function 

G{Ei,Ti;E2,T2) = a{Ei)a{E2) J dsids[ds2ds'2e-'^'^''^-''^-'^'^''^-'^^ 

y^VfAh - si)/a(ti - s[)f^{t2 - S2)fa{t2 " 4) ( ^2 (4) . ^^i (s'l )] [xa (sa) , a;i (si )] 

+A+[x2(4),a;i(si)]A+[xi(s;),X2(s2)]) , (52) 

refiects the correlation between the detection events in the two detectors. 

Eqs. fl5m52p apply to general paths in Minkowski spacetime. They simplify significantly, 
when we restrict to paths xiij) and X2{t), such that the Green's functions Air[xi(s), ^2(5')] and 
A[xi(s), X2(s')] are functions of s — s' only. Denoting the values of the Green's functions values as 
Air(s — s') and A"'"(s — s') respectively, Eq. (!52ll becomes 

G(Ei, n; E2, T2) = a{Ei)a{E2) J dSidS2Uri - 5i)/.(r2 - ^2) 

da;d2/e-'(^^+^^)"-'(^^-^^)^/'(7,(y2a;)(7.(y/\/2) 

X [A*p{S2 -Si + y)Ap{S2 ~Si-y) + A+(52 - 5i + x)A+{S, -S2 + x)] , (53) 

where 5*1 = (si + s[)/2, S2 = (s2 + S2)/2, x = ^{s[ — si + ^2 — S2) and y = s[ — Si — + S2, and 
we used the fact that ga{s'i — Si)ga-{s2 — S2) = ga{y/2x)ga{y/ V^)- 

We separate Eq. fl53|) into a sum of two terms, one involving the pair of Feynman propagators 
and one involving the pair of Wightman functions. Integration over x in the former term yields 
a proportionality factor e"'^ {^1+^2) . hence, this term is strongly suppressed. Integration over y 
in the latter term yields a proportionality factor a/ levrcr^ exp [— {E2 — -Ei)^cr^]. Since Eia » 1 
and £^20" >> 1, this factor is well approximated by 47r6{Ei — E2). This implies that non-trivial 
correlations exist only for the detection of particles with the same energy. Eq. (!53l) then becomes 

G{Eu n; E2, T2) = Ana{E,)a{E2)S{E, - E2) j dSidS2Mri - S,)Ut2 - S2) 

X J dxe-'(^i+^^)"^,(y2a;)A+(x + ^2 - ^i)A+(a; + ^1 - ^2). (54) 

3.2.2 Uniformly accelerated detectors 

Next, we consider two uniformly accelerated detectors separated by constant proper distance d in 
a direction normal to the direction of motion, that is, Xi(r) = (sinh(ar)/a, (cosh(ar) — l)/a, 0,0) 
and X2{t') = (sinh(ar'), (cosh(ar') — l)/a,o?, 0). The detector's clocks are synchronized so that 
xi(0) = (0, 0, 0, 0) and ^2(0) = (0, 0, d, 0). In this case, 

A [x,iT),X2ir')] = - sinh[a(r - r' - li - r) /2] sinh[a(r - r' - li + r) /2] ' ^^^^ 

where 

2 , fad\ , , 

r = — sinh — (56) 
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is the proper time that it takes a hght-ray emitted from a point of the path Xi(-) to reach a point 
of the path a;2(-)- 

We substitute Eq. (155|) into Eq. (15^ . employ Eq. (fTTj) and integrate over Si + S2, to obtain 



G{Ei,n;E2,T2)=a{Ei)a{E2)-^5{Ei-E2) J dS UAr ~ S)H{S), 
where 5 = 6*2 — Si, At = T2 — ti, and 

g-i(Ei+E2)x^^(y2x) 



(57) 



H{S) 



dx 



sinh[a(x — u)/2] sinh[a(x + u)/2] sinh[a(x + v)/2] sinh[a(x — v)/2] 



(5^ 



In Eq. dSSD, we wrote u = S + r and v = S — r. 

The function integrated in Eq. (!58|) has 4 poles, at x = ±m and x = ±v. The function 
strongly suppresses the integral Eq. (158!) . unless 1^1 is of order cr or smaller, or |f | is of order a or 
smaller. 

We find analytic expressions that approximate -^(5*) well in two regimes. The first regime 
corresponds to \u — v\ = 2r » a. The distance between the two detectors is sufficiently large, so 
that the time delay of signal propagation between them is much larger than the temporal resolution 
a. Then, the peaks at = u and |x| = f do not overlap. Hence, the integral Eq. fj58|) is well 
approximated by a sum of two terms, one corresponding to |x| = u and one to |a;| = v. To leading 
order in l/(cra). 



HiS) 



9a{V2' 



u 



sinh(aS') sinh(ar) 

gAV2v) 



-co—ie 
00— ie 



J sinh[a(x — m)/2] sinh[a(x + u)/2] 

^j.^-i(Ei+E2)x 



sinh(aS') sinh(ar 
Carrying out the integrations in Eq. (15^ . we obtain 

87r/a 



J sinh[a(x — v)/2] sinh[a(a; + v)/2] 



(59) 



H{S) 



(e^(^i+^2) _ X) sinh(aS') sinh(ar) 



g,{V2{S + T)) 



sin[(Ei + E2)(5 + r)] 



UV2{S - r)) 



sinh[a(S' + r)] 
sm[{Ei + E2){S-r)] 
sinh[a(S' — r)] 



(60) 



We substitute Eq. (160|1 into Eq. (157|) and integrate over S. In the Appendix, we evaluate this 
integral to leading order in [{Ei + E2)aY^ and (acr)~^. The result is 



^(^1, Ti- E2, T2) = -a{Ei)a{E2)8{Ei - E2) 



l^tanh 



■K{E-i+E2) 
2a 



e^(^^+^2) _ 1) sinh^(ar) 
x[^(Ar-r) + ^(Ar + r)] 



(61) 



The correlations are strongly peaked around At = ±r, i.e., information about the recording of a 
particle in a detector is transmitted with the speed of light and affects the response of the other 
detector. Since Eq. fl6Tl) holds for r >> a, it follows that ar » 1, and thus. 



-2ar 



tanh 



G{Ei, n; E2, T2) = -—a{Ei)a{E2)6{Ei - E2 
Zn 



■k{Ei+E2) 
2a 



e^{E,+E2) _ I 
x[U{AT-r) + UAT + r)] 



(62) 
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We observe that the correlations decrease exponentially with the time-delay r. 

The other physically interesting regime for which an analytic expression for the correlations 
is obtained, corresponds to small separation between the detectors, r « . In this case, the 
effects of retarded propagation are insignificant, and the trajectories followed by the two detectors 
are indistinguishable. Hence, in this regime, the two detectors may be viewed as independent 
subcomponents of a single macroscopic detector. At the limit r — )■ 0, Eq. fl58l) becomes 

^ ' y_oo-i.-sinh2[a(a;-5)/2]sinh2[a(x + 5)/2]' ^ ' 

The presence of the function g„ implies that the integral is strongly suppressed, unless S is of order 
a or smaller. To leading order in l/((Ta), we obtain 

^^^^1 ^ 167r^<,(V2S) 



a2(ef {£;i+i?2) _ i)sinh2(a^) 
X [{El + E2) cos[(Ei + ^2)^] - a coth(a^) sin[(Ei + ^2)^]] . (64) 

We substitute Eq. flM|) into Eq. and integrate over S. To leading order in [{Ei + E2)a]~^ 

and (a(j)~\ we obtain (see the Appendix) 



[El + E2 



\2 coth 



{Ei+E2)n 
2a 



GiEu n; E2, T2) = \ a{Ei)a{E2)S{Ei - E2) 2., ' ^ , U^r). (65) 



Second-order coherence functions. Eqs. ( 162|) and ( 165|) provide the two-time correlations for 
particle detection for a pair of uniformly accelerated detectors in the regimes ra » 1 and ra << 1 
respectively. It is convenient to express these correlations in terms of the second-order coherence 
function g^'^\ which is standardly employed in quantum optics. For a stationary field configuration, 
the second-order coherence function is defined as 

If g^'^^Ar) < g^'^\0), there is a higher probability for a simultaneous detection of a pair of 
particle, whence the particles are said to be bunching. If g^'^^Ar) > g^'^\0), the particles are said 
to be anti-bunching. When considering separated detectors, we have to take into account the effect 
of retarded propagation, so the natural definition of bunching is g^'^^Ar) < g^'^\r) for |Ar| > r. 

Also, the statistics of particle detection are sub-Poissonian if ^'^^•'(O) < 1, Poissonian if ^'^^•'(O) = 
1, and super- Poissonian if ^''^^•'(O) < 1. 

The intensity / and the intensity correlations are defined by Eq. and Eq. fl33p , respectively. 
Using Eqs. flSTj) . f l62|) and f l65|) . we obtain 



X 2 2arfdEE''a(Eye-'^sech\7iE/a) 

J dEa{E)Ey{e'-^ - 1)1 
X [faiAr - r) + /^(Ar + r)] , for r» a (67) 

N f dEE'^a(E)'^e~~^cosech.'^(7TE/a) „ , , , , 

g^^\AT) = 1-^^-^ — ^2^/,(Ar),for r = 0. (68) 

fdEa{E)Ey{e'-^ - 1) 
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An important special case is that of a macroscopic two-level detector, namely, a detector 
recording only a narrow frequency band [E — ^AE,E + ^Ai?], with << AE « E. In this 
case, 

g^^\AT) = l--^^e-2-tanh(7rE/a)[^(Ar-r) + ^(Ar + r)],for r »a (69) 
<?(2)(Ar) = l-||/^(Ar),for r = 0. (70) 

It is evident from the equations above that particles recorded by a uniformly accelerated de- 
tector anti-bunch and satisfy sub-Poissonian statistics, at both regimes ra » 1 and r — )■ 0. 

3.2.3 Static detectors in a thermal bath 

In order to examine whether the relation between uniform acceleration and temperature holds, 
we must construct the second-order coherence function corresponding to a pair of static Unruh- 
DeWitt detectors when the field is in a thermal state of temperature T. To this end, we compare 
the positive frequency Wightman function for a thermal state Eq. (jH]) with the one corresponding 
to a uniformly accelerated observer Eq. fl55l) . One obvious difference is the presence of the delta 
functions in Eq. ( H^ : these, however, do not contribute to Eq. ( 15^ because they vanish along 
the real axis. We denote the regular part of the thermal Wightman function Eq. (jH]) as A^ and 
the Wightman function for the accelerated observers Eq. (!55l) . evaluated for a = 27r//3 as Aa=27r//3- 
Then, 

/3 f2nr\ 

= ^ Si^ll — J ^a=2n/p. (71) 

In Eq. (I7T]) . r refers to the time delay due to retarded propagation. For static detectors in a 
thermal bath, r coincides with their spatial separation. Thus, the difference in the two Wightman's 

function is a multiplicative factor e(r) = ^ sinh j that tends to unity at the limit r — )■ 0. 

This implies that the function G{Ei,ti] E2,T2) for static detectors in a thermal field is obtained 
from Eqs. (162|) and (!65|) by setting a = 27r//3 and multiplying by a factor e(r)^. 
In the regime r » a, 

G{Ei,n;E2,T2) = ~-^a{Ei)a{E2)5{Ei - E2) 

X ^5^^|g^[/.(A.-.).«A...)l. (72) 

and for r = 0, 

r(T? p ^ (-^1 + E2? ^ , coth[/3(Ei + E2)/4] ^ 

G(Ei, n; E2, T2) = — a{Ei)a{E2)S{Ei - E2) ^^^e.+e^) _ i /-(^^)- (^3) 

Then, the second-order coherence functions are 

^ vr jdEEMEfe-^''^ech\(5El2) 

^ ^ 2r2 [j dEE^a{E)/{eP'^ 

X [/^(Ar-r) + /^(Ar + r)],forr »(7 (74) 

f2WA N . 'RjdEE^a{Efe-P'^cosech'{PE/2)^^^ 
- 2[/.eU(W^-1)P 
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For a macroscopic two-level detector, 



^7(2)(Ar) = 1- -^^^^ tanli(/3E/2) [^(Ar-r) + ^(Ar + r)], for r»a (76) 
^(2)(Ar) = l-||/,(Ar),forr = (77) 

Comparing Eq. fl75|) and Eq. fl68|) we see that for r = 0, the intensity fluctuations recorded by 
uniformly accelerated detectors and the ones recorded by static detectors in a thermal bath are 
fully equivalent, given Unruh's relation between temperature and acceleration T = However, 
this equivalence is lost for separated detectors, in the regime r » a. While in both cases we have 
particle anti-bunching, the dependence of the correlations on the time delay r differs significantly. 
A pair of inertial detectors in a thermal bath has long range correlations (that fall-off with r~^), 
while a pair of accelerated observers has short-range correlations that drop exponentially with r. 

It is important to emphasize that the difference above arises at the level of the field's two-point 
functions, which are expected to enter into the correlations for any observable defined with respect 
to the two separated detectors. One therefore expects that the correlations of separated accelerated 
observers will generically deviate from the thermal ones. 



3.3 Relation to the Glauber detector 

The fact that the second-order coherence function ( 1751) for the thermal state exhibits anti-bunching 
(g^'^^Ar) > g^'^\0) for r = 0) might be viewed as contradictory to the well-known fact that thermal 
photons bunch [T7]. The difference here is that we consider Unruh-DeWitt detectors, while the 
bunching behavior of photons in quantum optics is derived for Glauber detectors. 

In Glauber-type detectors, the n-time detection intensity correlation function at different space- 
time points xi,X2, . . ■ ,Xn is proportional to 

G(")(xi, 0:2, ...,Xn) = Tr f0(~)(xi) . . . 4>^-\xn)^^+\xn) . . . ^^'\xi)po) , (78) 



where 0*^"-' and 0^^-' are respectively the positive- and negative frequency components of the 
Heisenberg-picture scalar field 0(x), and po is the field's density matrix. 

The correlation functions of the Glauber detectors can be obtained from the general theory of 
Sec. 2, in particular, Eq. The interaction Hamiltonian for a Glauber detector is of the form 

Hj = cr+0(-)(a;) + a-0W(x), (79) 

where a^Po = Pqo^ = 0. 

It is possible to bring the Unruh-Dewitt interaction Hamiltonian into the Glauber form through 
a redefinition of variables and the consideration of 'dressed' states for field and detector [H]. For 
this reason, we believe that the interaction Hamiltonian does not constitute the most important 
difference between the two types of detector. As explained below, the most important difference 
is the scale of the temporal coarse-graining parameter a. 

The two-time probability distribution Eq. ( 132|) for a pair of detectors then becomes 



p{Ei,Ti] E2,T2) = a{Ei)a{E2) J dsids[ds2ds'2 x e 



■iEi{s[-si)-iE2{s'2-S2) 



X V/<7i(tl - Si)/^i(ti -S'i)/^2(t2 - 52)7^2(^2 - 4) 



xTr 



PoA (0(-)[xi(.i)]0^-^k2(4)]) r [^^^\x2is20^\x,{s^)])] , (80) 
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where Xi(-),X2(-) are the detector's world-hnes. 

Eq. f lSU]) leads to a two-time probabihty proportional to C^") of Eq. fl75]) at the limit cr — )■ 0. 
This means that the resolution timescale a of the detector must be smaller than any time-scale 
appearing in the correlation function 



Tr 



(81) 



In contrast, an Unruh-Dewitt detector recovers the Planckian spectrum only if aa » 1, i.e., if the 
resolution timescale is much larger than the characteristic timescales of the correlation function. 



4 Conclusions 

In this paper, we constructed new models for macroscopic particle detectors, where the time of 
transition is treated as a quantum observable. These models allow for a precise treatment of 
temporal fluctuations and correlations in particle detection. For Unruh-DeWitt detectors, moving 
along generic trajectories in Minkowski spacetime and coupled to a quantum field, we obtained 
a causal, and approximately local-in-time expression for the detector's response: the response 
function depends only on properties of the path that are defined at time t, within an accuracy 
allowed by the detector's resolution. This implies that Unruh's relation between acceleration and 
temperature is quite robust and extends beyond the case of uniform acceleration: slow changes 
of the detector's acceleration result to corresponding changes in the temperature of the recorded 
Planckian spectrum. 

Our most important result is the explicit calculation of the two-time probabilities and of the 
intensity fluctuations for uniformly accelerated detectors. We find that that intensity fluctuations 
recorded by a single detector are thermal: the second-order coherence function of a uniformly 
accelerated Unruh-DeWitt detector is identical to that of a static detector in a thermal bath. This 
implies that the records of particle detection by a single Unruh-Dewitt detector cannot distinguish 
whether the detector is accelerated, or whether it is static in the presence of a thermal bath. In 
contrast, the intensity correlations for a pair of Unruh-DeWitt detectors, with the same linear 
acceleration, but moving along different trajectories are non-thermal. Hence, records of particle 
detection from a pair of detectors, at rest with respect to each other, do distinguish whether the 
detectors are in an accelerated reference frame or whether they are static and in contact with the 
thermal bath. This result greatly strengthens the view of the Unruh effect as a fundamentally 
local phenomenon [5^ . 

It is important to note that the existence of non-thermal correlations cannot be inferred from 
the quantum-field theoretic treatment of the Unruh effect (Fulling- Rindler quantization). In that 
treatment, the restriction of the Minkowski vacuum in one Rindler wedge is a thermal state with 
respect to the Rindler time coordinate that is associated to accelerated observers. If this restricted 
state were interpreted as the quantum state of the field in the accelerated reference frame, then one 
would expect that all field observables with support on one Rindler wedge would exhibit thermal 
properties in the accelerated frame, irrespective of whether these observables are local or not. Our 
results suggest that the Fulling-Rindler quantization should not be taken as a literal construction 
of the quantum state in an accelerated reference frame [IQ]. Indeed, the interpretation of the 
Fulling-Rindler vacuum has no bearing on the, fundamentally local, physics of the Unruh effect 
|39] . In contrast, the detectors employed here define physically meaningful observables for the 
quantum field, with a natural interpretation in terms of accelerated frames of reference. Hence, 
the conclusion that the particle-detection correlations are not thermal provides novel informa- 
tion that cannot be extracted from the procedure of Fulling-Rindler quantization, at least in a 
straightforward way. 
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The question then arises how one should describe the quantum field state in the reference 
frame of accelerated observers. In our opinion, the most promising approach would involve the 
framework of local quantum physics |42J, with a construction of an algebra of local observables 
and the associated space of states. In order to explain this point, we note that the n-time joint 
probabilities Eq. can be expressed in terms of positive operators P„(i?i, ti, E2, ^2, • • • En, tn) 
on the Hilbert space of the quantum field, by tracing out the detectors' degrees of freedom. When 
considering detectors moving along the same trajectory, the operators P„(i?i, ti, £'2, ^2, • • • -^n, in) 
define a subset of local observables with thermal expectation values. This description fits naturally 
with the idea of locally thermal states, defined in Ref . [13] . In this work, states of local thermal 
equilibrium are defined if the expectation value of a sufficiently large number of thermal observables 
in the vicinity of a spacetime point x coincide with the expectation value of a globally KMS state 
at this point. This analogy suggests the possibility of describing the quantum field state in an 
accelerating reference frame as a local equilibrium state, where the local temperature at each point 
is proportional to the proper acceleration in this frame, but with non-local correlations between 
different spacetime points. 

The fact that there may exist significant information about the correlation functions of the 
field that is not captured by transformations between different field vacua is an important obser- 
vation that could be relevant to the information-loss paradox in black holes. There are, of course, 
significant differences between the physical set-ups of the Unruh and the Hawking effects, so the 
results of the present paper cannot be directly applied to black holes. Nonetheless, the point of 
principle is intriguing. The information-loss paradox arises from the assumption that the final 
state of the field at late times is globally Gibbsian, hence, mixed. If, however, the state is non 
Gibbsian, but only locally thermal, it is at least conceivable that the information of the initial 
state may persist throughout the process of black hole formation and evaporation, in the form of 
non-local correlations of the field variables. 



A Evaluation of integrals 

Here, we provide the derivation of Eqs. (16 ip and (16 5 p from the evaluation of integrals of the form 
Eq. dSZl). 

In order to derive Eq. f l6T]) . we must compute the integral 



dSUS-/\T)H{S), (82) 
where H{S) is given by Eq. ( pUj) . This means that we must evaluate integrals of the form 

J smh(a5) smh[a(D ± r)J 

For concreteness we compute the integral above with the (-) sign. The integrand is peaked around 
S = At and S = r. We change the integration variable to x = S* — r, so that 



J smh[a(x + r)J smh(ax) 



where v = At — r. We set A{x) = fa{x — v)ga-{V2x)/ sinh[a(x + r)] and B{x) = Denoting 
by A{k) and B{k) the Fourier transforms of A{x) and B{x) respectively, Eq. (|8^ becomes 

-A{k)B*{k). (85) 
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For r >> 0", the smh[a(a; + r)] term is well approximated by sinh(ar), and, consequently, A{k) is 
proportional to e"*^ This means that values of >> cr^^ are suppressed in the integration 
Eq. (!85|) . The function B{x) is oscillatory with a frequency of order E; hence, to leading order in 
(Ea)-^, we can substitute IS*{k) with 5(0) in Eq. §5^. Thus, 



We compute 



Hence, 



dk 

Ji ~ 5(0) / —A{k) = B{0)A{0). 
Zn 



(86) 



sm(Ex) TT 
rfa; — — — - = — tanh 
sinh(aa;) a 



7tE 

'2a 



17) 



Ji 



TT 



asinh(ar) 



tanh 



ttE' 



(88) 



from which Eq. fl6Tl) follows. 

The derivation of Eq. ( l65l) proceeds along similar lines. It involves an integral of the form 



J2 



dS 



sinh {aS) 



f^{S - AT)g^{V2S) [Ecos{ES) - acoth{aS) sm{ES)] 



(89) 



Using the same arguments as before, we find that to leading order in {Ea) ^, J2 = /o-(^t)-B(0), 
where now 



dS 



E cos{ES) - a coth(a^) sm{ES) _ ttE^ coth (|f ) 



sinh^(aS') 



2a2 



(90) 



Hence, 



7lE^COth(^) , ^ 

.h = ^/.(A-) 



(91) 



from which Eq. f l65|) follows. 
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